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Abstract — The paper examined the cyclical trends, seasonal variation and seasonal adjustment of global 
CO 2 emission from 1970 to 2018 through the application of Hamilton regression filter model. ARIMA (4,0,0) 
forecasting model for 2030 has been added with the Hamilton filter model and observed that the new model 
is stable, stationary and significant in which volatility is being minimised and the heteroscedasticity problem 
is totally disappeared. 
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I. INTRODUCTION 

A series contains both the cyclical and seasonal behaviours. 
A seasonal behaviour is very strictly regular, meaning there 
is a precise amount of time between the peaks and troughs 
of the data. Cyclical behaviour on the other hand can drift 
over time because the time between periods isn't precise. A 
series with strong seasonality will show clear peaks in the 
partial auto-correlation function as well as the auto¬ 
correlation function, whereas a cyclical series will only 
have the strong peaks in the auto-correlation function. 
Patterns of peaks and troughs occur in both cyclical and 
seasonal analysis wherein seasonal trends the period 
between successive peaks (or troughs) is constant. But in 
cyclical trends the distance between successive peaks is 
irregular. 

The decomposition of macroeconomic time series with 
trend and cyclical components is crucial to many 
macroeconomic concepts such as potential output, natural 
interest rate, share prices and inflation rate and so on. 
Econometric analysis requires filtering techniques which 
cater data sequences having short and strong trends. Eilter 
can be used to decompose a time series into trend, seasonal 
and irregular components. 

H.P.Eilter (1997) has been used as a low pass smoothing 
filter in a numerous macroeconomic investigation setting 
smoothing parameter to certain arbitrary conventional 
values. It can decompose trend and cycles of the time series 
data to produce new time series such as potential GDP and 
output gap that are useful in macroeconomic discussion and 
in some recent public debate. It also used ARM A 


forecasting technique that deliver the required forecasted 
values that are needed in H.P. algorithm. The cyclical 
component obtained by H.P.Eilter is unpredictable which 
indicates that we cannot do any investigation into the 
cyclical component because they are independent. 

Many business cycles analysis, volatility of stock price 
indices and unemployment trends and cyclical patterns have 
been being studied through Baxter and King (1999) and 
Christiano and Eitzgerald (2003) models of decomposition 
of trends and cycles through filter till now. 

Therefore, the author endeavours to examine the patterns of 
trends and seasonal and cyclical fluctuations of global CO 2 
emission during 1970-2018 by applying the methodology of 
Hamilton Eilter model(2018).The paper also finds the 
projected rise of global CO 2 emission in 2050 and examined 
the forecasting ARIMA(p,d,q) model incorporating in the 
Hamilton regression filter model up to 2030. 

II. SOME IMPORTANT STUDIES 

Bhowmik(2018) studied that global CO 2 emission has been 
rising at the rate of 2.19% per year and global per capita 
CO 2 emission has been catapulting at the rate of 0.58% per 
annum during 1960-2015.Global CO 2 emission contains 
four upward structural breaks in 1968,1976,1988 and 2004 
and per capita emission has two upward structural breaks in 
1969 and 2004 respectively. World CO 2 emission is 
positively related with global GDP and GDP per capita 
significantly during the specified period.There is no 
cointegration between world GDP and world CO 2 emission 
during 1960-2015. 
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Hamilton filter (2018) claims that H.P.Filter (1997) creates 
spurious dynamic relation, fails to eliminate spurious 
predictability,produces values for the smoothing parameter 
vastly at odds with common practice. It is not also a sensible 
approach for a random walk. Therefore, author uses 
regression variable at date t+h on the four most recent 
values as of date t which is a better alternative to detrend 
and to extract cyclical components. Hamilton also assumes 
random walk, then cyclical component is estimated as yt+h - 
yt which is predictable and perform much better than the 
H.P.Filter model. 

Ciais et al.(2018) examined that since the beginning of the 
Industrial Era, fossil fuel extraction from geological 
reservoirs, and their combustion, has resulted in the transfer 
of significant amount of fossil carbon from the slow domain 
into the fast domain, thus causing an unprecedented, major 
human-induced perturbation in the carbon cycle during 
2000-2009 which was observed a major change over 
preindustrial period. Emission of carbon from fossil fuel 
reserves, and additionally from land use change is now 
rapidly increasing atmospheric CO 2 content. 

Schuler (2018) applied Hamilton filter in the credit to GDP 
gap in USA taking data of GDP from Ered during 1947Q1- 
2017Q1 and data of credit GDP ratio from BIS during 
1952Q1-2016Q4 which produced more robust cycle 
estimates than H.P. filter. 

Bhowmik(2020) examined the growth rate of GDP in India 
during 2011Q4 -2019Q4 using quarterly data in explaining 
the trend,seasonal patterns and cycles in context of 
Hamilton regression filter model.The paper concludes that 
the TRAMO/SEATS of Hamilton filter produce clear 
seasonal variation showing v shaped volatility and cycles 
confirmed two peaks. ACF and PACF of Hamilton residual 
series verified the nature of seasonality. Residuals suffer 
from heteroscedasticity. ARIMA( 1,0,0) process of 
Hamilton filter for forecasting 2030 converges to stable and 
stationary processes. 

Even,Hodrick(2020) applied ARIMA model in Hamilton 
regression filter in decomposing cyclical components of US 
GDP during 1947Q1-2019Q1 using simulation technique 
and found that coefficient of AR( 1) is greater than one and 
coefficients of AR(2) and MA(l)and MA(2) are less than 
one and the estimated equation is significant. 

Oliver and Peter(2020) stated that global CO 2 emission has 
increased by 1.6% per year during 1970-2003 which 
accelerated to 3.2% per year during 2003-2011 using the 
data set of EDGAR V4.2ET2010. 


EDGAR(2020) stated that 2015 was the hottest year among 
the 16* warmest year during 1998-2015 but there was a 
structural break since annual growth rate of CO 2 from 2012- 
2015 decreased especially in China, USA, EU, Canada and 
there was increased decoupling of CO 2 emission and GDP. 

Lindsey(2020) reported that during 1750-2019,the 
atmospheric concentration of CO 2 in ppm has been steadily 
increasing but the CO 2 emission has been catapulting along 
with small cycles. Author mentioned that CO 2 concentration 
was never exceeded over 300ppm until preindustrial 
revolution, but in 1958 it reached to 315 and in 2015 it 
accelerated to 400ppm.If global energy demand continues 
to grow,the CO 2 emission is projected to exceed 900ppm by 
this century(2100). 

III. METHODOLOGY AND SOURCE OF DATA 

Semi-log regression model is applied to find linear trend 
which is given by, 

Log(xi)=aH-btH-Ui 

Where Xi= variable to be estimated, a and b are constant,t 
=time,Ui=random error and Ul,2,.n. 

Then Cochrane and Orcutt(1949) model was applied to find 
the projection of global CO 2 emission for 2030. 

Hamilton (2018) regression filter model has been calculated 
in the following manner to find the decomposition of 
cyclical trends, cycles and seasonal variations. 

yt+8=ao+aiyt+a2yt-i+a3yt-2+a4yt-3+Vt+8 
Or, Vt+8=yt+8+do+diyt+d2yt-i+d3yt-2+d4yt-3 
So, yt=ao+ aiyt-8+a2yt-9+a3yt-io+a4yt-ii+Vt 

Therefore,Vt= yt-(do+diyt-8+d2yt-9+d3yt-io+d4yt-ii) where di 
are estimated and y is the variable to be estimated. 

Vt+h=yt+h -yt is the difference i.e. how the series changes over 
h periods. For h=8, the filter 1-L^ wipes out any cycle with 
frequencies exactly one year and thus taking out both long 
run trend as well as any strictly seasonal components. 

It also applies random walk: yt=yt-i-£t where d=l and 
C0t^=Et+h+Et+h-1+.+Et+1 

Regression filter reduces to a difference filter when applied 
to a random walk. Hamilton suggested h=8 for business 
cycles and h=20 for studies in financial cycles. Regression 
Vt converges in large samples to ai=l and all other aj 
=0.Thus, the forecast error is Vt+h=yt+h -yt. 

The equation Vt can be decomposed into trend, cycle and 
seasonally adjusted through SEATS/TRAMO or STL or 
census X-13 packages. The Box and Jenkins(1976) model 
The forecasting of ARIMA(p,d,q) has been applied in the 
Hamilton regression filter model for prediction in 2030. 
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of ARIMA(p,d,q) can be estimated as below. 

Vt=a+PoiVt-i+Et+PiEt-i+et 

where Vt=Hamilton regression filter residual, a is constant, 
Poi andPiare the coefficients of AR and MA. AR terms and 

MA terms are Vt-iand Et-irespectively where i=l,2.n 

,andet= random error/residuals. For convergence criteria, 
Poi andPishould be less than one and the values of AR and 
MA roots must be less than one for stability and stationary 
condition. 

Data of global CO 2 emission in Giga Ton from 1970 to 
2018 were collected from the World Bank. 


IV. OBSERVATIONS FROM THE MODELS 

Global CO 2 emission has been catapulting at the rate of 
1.74% per annum significantly during 1970-2018 which is 
estimated by the trend line given below. 

Log(xi)=2.74lH-0.0174tH-u, 

(43.94)*(9.90)* 

R2=0.99,DW= 1.31, AIC=-235.31 ,SC=-231 .56,xi=global 
CO 2 emission in giga ton,*=significant at 5% 
level.,Ui=residual. 

In 2018,the global CO 2 emission was 36.10gigaton which 
is projected to rise 63.503Gigaton in 2050 which is 
depicted in Figure 1 below. 



year 


Fig.l: Global CO 2 emission projection 


Source: Plotted by author 

Actually, the global CO 2 emission fulfils the quadratic 
function properties which indicated that the trend line is 
non-linear and the estimated equation is given below. 

Log(xi)=2.748H-0.013tH-0.000105t2H-Ui 
(159.11)* (8.31)* (3.41)* 

R2=0.98,F=1129.95*,DW=0.314,*=significant at 5% 
level, Ui=residual. 


This non-linear trend is significant at 5% level having 
autocorrelation problem. 

The global CO 2 emission has the property of having five 
structural breaks in its path from 1970 to 2018 in 
1977,1987,1996,2003 and 2010 all of which are significant 
at 5% level and all showed upward rising. It is depicted in 
Figure 2 below. 
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Source: Plotted by author 


- Residual - Actual - Fitted 


Fig.2: Structural breaks of global CO 2 emission 


Hamilton regression filter model of global CO 2 emission 
during 1970-2018 is estimated below where the adjusted 
sample has fallen during 1981-2018 where total 
observations are 38 after adjustment. 

Log(xi)t=-0.2489-F0.69541og(xi)t-8-F0.28721og(xi)t-9 
(-1.59) (1.50) (0.38) 


respectively and all the coefficients are less than one.The 
residual Vt is calculated as 

v,=log(xi) -[-0.2489-F0.69541og(xi)t-8-F0.28721og(xi)t-9- 
0.57411og(xi),.io+0.72631og(xi),_ii] 

This equation can be treated as Hamilton filter model of 
global CO 2 emission. 


-0.57411og(xi),.io-F0.72631og(xi)t-ii-FVt 
(-0.72) (1.48) 

R2=0.94,F=129.50*,DW=0.283,AIC=-2.71,SC=-2.49 

The Hamilton regression filter is found as a good fit with 
high and significant F values with low AIC and SC 


The estimated Hamilton regression filter equation is plotted 
in Figure 3 where fitted and actual lines are shown as 
upward rising and the residual consists of 7 peaks and 
troughs respectively whose cyclical trend is diverging away 
from equilibrium with rising volatility. 



Source-Plotted by author 
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Hamilton filter residual of global CO 2 emission during 
1970-2018 has been decomposed into trend, cycle, seasonal 
variation, remainder, seasonally adjusted cycle respectively 
in a panel of diagrams in Figure 4 below. In panell,the 
Hamilton filter residual is shown as a cyclical pattern 
having seven peaks and troughs. In panel 2,the smooth 


cyclical trend is plotted which minimises cycles showing 
two peaks and troughs where amplitude is widening. In 
panel 3,the seasonal variation is shown where its patterns 
are v shaped. In panel 4,the remainder is presented which is 
fluctuating and in panel 5,the seasonally adjusted cycle is 
visible which looks as like as residual line. 
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Fig.4: Decomposition of Hamilton regression filter of global CO 2 

Source-Plotted by author 


In the Figure 5, the Hamilton filter of global CO 2 emission, 
its cycle, trend and seasonal adjustment is plotted in the 


composite figure. The Hamilton filter trend line exemplified 
with two peaks and troughs which are diverging with bigger 
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amplitude. There is no wide difference between seasonally 
adjusted and the actual residual of the filter of global CO 2 , 


although both of which are cyclical with upswings and 
downswing phases. 



Source-Plotted by author 


The seasonal fluctuation can be confirmed by the 
application of finding auto correlation functions and partial 
autocorrelation functions of the Hamilton regression filter 
residual series of the global CO 2 emission from 1970 to 
2018.The ACF have been declining up to lag 4 and then it 
reached negative which are increasing up to lag Hand 
started to diminish up to lag 14 and then it moved up 
towards positive values from lag 15.Similarly PACF 


declined up to lag 2 and then turned to negative up to lag 4 
and from lag 6 the negative values increased up to lag 10 
and reached positive values up to lag 13 and then fell 
negative. The patterns of ACF and PACF indicate that the 
seasonal fluctuation of global CO 2 emission is confirmed 
since the probabilities of Q stat are significant. In Figure 6,it 
proves clearly. 


Autocorrelation 

Partial Correlation 

AC PAC 

Q-Stat 

Prob* 

1 


1 


1 0 688 0 688 

19 435 

0 000 

1 


1 


1 

2 0.483 0.018 

29.271 

0.000 

1 


1 


1 

3 0.337 -0.002 

34.211 

0.000 

1 

1 1 



1 

4 0.037 -0.378 

34.271 

0.000 

1 1 

1 

1 


1 

5 -0.095 0.020 

34.686 

0.000 

1 ■ 

1 

1 1 


1 

6 -0.175 -0.047 

36.148 

0.000 

!■ 

1 

1 1 


1 

7 -0.292 -0.094 

40.328 

0.000 


1 

1 1 


1 

8 -0.304 -0.069 

45.015 

0.000 


1 

1 1 


1 

9 -0.306 -0.075 

49.917 

0.000 


1 

1 ■ 


1 

10 -0.373 -0.200 

57.465 

0.000 


1 

1 1 


1 

11 -0.378 -0.136 

65.516 

0.000 

!■ 

1 

1 

1 

1 

12 -0.284 0.088 

70.219 

0.000 

1 ■ 

1 

1 

1 

1 

13 -0.173 0.112 

72.047 

0.000 

1 1 

1 

1 


1 

14 -0.059 -0.015 

72.268 

0.000 

1 

1 1 

1 1 


1 

15 0.066 -0.036 

72.552 

0.000 

1 

1 1 

1 ■ 


1 

16 0.076 -0.182 

72.954 

0.000 


Fig.6: AC and PAC of Hamilton filter residual 
Source-Plotted and calculated by author 

Among the best 20 models of automatic forecasting of regression filter models of global CO 2 emission,the 

Autoregressive Integrated Moving Average of Hamilton ARIMA(4,0,0) model is selected for forecasting for 

ISSN: 2456-1878 

https://dx.doi. 0 rg/lO. 22 i 6 l/iieab. 55.9 


1255 
























InternationalJournal of Environment, Agriculture and Biotechnology, 5(5) 
Sep-Oct, 2020 / Available: https://iieab.com/ 


2030.By applying ARMA maximum likelihood method 
after adjusting sample during 1981-2018,the estimated 
equation is observed as below. 

v,=0.0022lH-0.8635vt-iH-0.0452vt-2H-0.4387v,-3 
(0.127) (4.98)* (0.18) (2.17)* 

-0.6406vt-4+0.000453o2t 
(-4.85)* (4.11)* 

R2=0.84,F=35.65*,DW=2.084,AIC=-4.44,SC=-4.18,AR 
roots=0.89±0.32i,MA roots=-0.46±0.70i 


coefficient of o^ is minimum and significant along with 
other three coefficients.Even,R^ is high,F is significant,AIC 
is minimum and the roots of AR and MA are less than one 
which revealed that the both AR and MA processes are 
convergent,stable and stationary. 

In Figure 7,the forecast AR1MA(4,0,0) is depicted neatly 
where the forecast value for 2030 is (0.044088) in 
comparison with the value of (-0.0950) at 2018 and forecast 
line is moving upward up to 2027 and then it is declining 
although its amplitudes were increasing up till 2010. 


The AR1MA(4,0,0) forecasting model of Hamilton residual 
is highly significant showing minimum volatility where the 



year 


Source- Plotted by author 


Fig. 7: ARIMA forecast in Hamilton filter model 


The heteroscedasticity test for residual of AR1MA(4,0,0) by 
applying ARCH(4) revealed that nR^=3.0166 because 
n=34,R^=0.088872 which is accepted by Chi-square(4) test 
whose probability is found as 0.55 ( where F=0.7058 whose 
probability is 0.59) after adjusting the sample during 1985- 
2018 which implied that no heteroscedasticity is accepted at 
null hypothesis. 

V. LIMITATIONS AND SCOPE FOR FUTURE 

STUDIES 

The paper did not compare the observations of 
decompositions of Hamilton regression filter model with 
the Hodrick-Prescott Filter model.If the quarterly or 
monthly data of CO 2 emission were taken then the 
observations from Hamilton filter will be more perfect and 
analytic. The details of autocorrelation and partial 


autocorrelation functions are not shown step by step of all 
lags. Application of Hamilton filter through ARJMA(p,d,q) 
forecast model may produce many debates. The 
clarifications of explaining amplitudes were not included in 
this analysis.Fven,the explanation of volatility through 
appropriate GARCH model is also absent 
here.Therefore,there is enough scope of future research in 
the context of this paper. 

VI. MODEL DEMANDS POLICIES 

Based on Business as Usual trend,global CO 2 emission is 
projected to increase some 43.08 billion metric ton in 2050 
in comparison to 35.3 billion metric ton in 2018.By 
2100,CO2 concentration will reach between 540-970ppm as 
compared to 368ppm in 2000.(https://www.acer- 
acre.ca/resources/climate-change-in-context/introduction- 
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2/global/scientific-projections/projections-for-carbon- 
dioxide).By 2100,cumulative global CO 2 emission will 
exceed 1 trillion tons of carbons (3.67 trillion of CO 2 ) 
threshold which according to IPCC will raise the earth’s 
surface temperature to 2°C above the pre-industrial 
minimum and trigger dangerous interference. 
(http://euanmearns.com/global-co2-emissions-forecast-to- 
2100 /). 

Paris Climate Agreement targets to limit temperature rise 
with in 1.5°C that might increase global GDP by 60-65% by 
2030.The agreement claims to reduce CO 2 emission by 20% 
along with 20% increase in renewable energy and energy 
efficiency respectively. As a whole, GHG emissions 
intensity must be reduced by 33-35% below 2005 level 
within 2030. So, to limit global temperature rise to 
1.5°C,emissions must be below 25gigatons by 2030.As on 
2020,February,194 states with EU agreed the Paris 
Agreement where 188 countries plus EU represent more 
than 87% of GHGs. China had committed to cut GHGs by 
26-28% by 2025 and USA assured to cut through Clear 
Power Plan. India is also committed to reduce the emissions 
intensity of its GDP by 33 to 35 percent by 2030 from 2005 
levels. India’s emission intensity has reduced by 21% over 
the period 2005-2014. By 2030, India’s emission intensity 
is projected to be even lower—in the range of 35 to 50 
percent, ( https://www.nrdc.org/experts/aniali-iaiswal/road- 
paris-indias-advancement-its-climate-pled) . 

VII. CONCLUSIONS 

The paper concludes that the global CO 2 emission has been 
catapulting at the rate of 1.74% per annum significantly 
during 1970-2018 and its projected value of 2050 is 
63.503gigaton in comparison to 36.10 gigaton in 
2018.Hamilton regression filter showed that global CO 2 is 
cyclical in patterns having 7 peaks and troughs but the 
smooth cyclical trend consists of two distinct peaks and 
troughs which are widening in the amplitude and the 
seasonal variation is v shaped when it is decomposed by 
STL method. The seasonal variation is confirmed by the 
patterns of ACE and PACE whose Q stat are significant at 
5% level. Automatic forecasting ARIMA(4,0,0) model for 
2030 is incorporated with Hamilton filter and the result is 
appeared as convergent, stable and stationary with 
minimum volatility. This model is free from 
heteroscedasticity problem. 
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